Identifying Deficiencies of Standard Accretion Disk Theory: 
Lessons from a Mean-Field Approach 

A. Hubbard 1 ' 2 , E.G. Blackman 1 ' 2 

1. Dept. of Physics and Astronomy, Univ. of Rochester, Rochester, NY 14627; 2. 
Laboratory for Laser Energetics Univ. of Rochester, Rochester NY, 14623 

ABSTRACT 

Turbulent viscosity is frequently used in accretion disk theory to replace the 
microphysical viscosity in order to accomodate the observational need for in- 
stabilities in disks that lead to enhanced transport. However, simply replacing 
the microphysical transport coefficient by a single turbulent transport coeffi- 
cient hides the fact that the procedure should formally arise as part of a closure 
in which the hydrodynamic or magnetohydrodynamic equations are averaged, 
and correlations of turbulent fluctuations are replaced by transport coefficients. 
Here we show how a mean field approach leads quite naturally two transport 
coefficients, not one, that govern mass and angular momentum transport. In 
particular, we highlight that the conventional approach suffers from a seemingly 
inconsistent neglect of turbulent diffusion in the surface density equation. We 
constrain these new transport coefficients for specific cases of inward, outward, 
and zero net mass transport. In addition, we find that one of the new transport 
terms can lead to oscillations in the mean surface density which then requires 
a constant or small inverse Rossby number for disks to maintain a monotonic 
power-law surface density. 

Subject headings: accretion discs-planetary systems:protoplanetary discs-turbulence 



1. Introduction 



Accretion disks are ubiquitous around stars and compact objects (jFrank et al.l 120021 ) . 
Those around black holes produce some of the most luminous objects in the universe, while 
those around young stars provide the material from which planets form. Thin astrophysical 
accretion disks are primarily angular momentum supported and thus to be accreting, require 
a mechanism to export angular momentum. Observations suggest that micro-physical vis- 
cosity acting on the differential rotation in a near-Keplerian disk is too weak to accommodate 
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needed time scales and luminosities, given mass density constraints (jPringld Il98ll ) and so 
attention has focused on turbulence to enhance the needed transport. While outflows are 
also likely important, the focus on turbulent transport raises three issues: (1) What are the 
minimum properties that turbulence must have to transport angular momentum outward? 
(2) What causes the turbulence? (3) What practical set of equations captures the correct 
physics? 

Turbulence was first inc orporated into a practical analytic formalism by Shakura-Sunyaev 
(IShakura fc Syunyaevl Il973l ). There, the only role of turbulence emerges as a simple re- 
placement of the micro-physical viscosity in the momentum equation with an enhanced 
turbulent viscosity. One way that this immediately falls short is that the current leading 
candidate for turbulent transport is the magneto-rotational instability (MRI) which involves 
magnetic fields not considered in Shakura-Sunyaev formalism. That being said, the MRI 
drives turbulence and one might hope that the resulting transport could be modeled via a 



Shakura-Sunyaev formal i sm. Evidence is emerg ing that this is not the case (IPessah et al. 



20071 . IPessah et al.l 12008k iBlackman et al.l 120081 ) and there presently remains a disconnect 
between the insights gained from numerical simulations and a practical formalism that can 
both accommodate the M!RI ph ysics and be of value to phenomenological modelers (e.g. 
Pessah et al.ll2007LIOgilviell2003h . 



As part of a long term effort to develop a formalism which captures the minimum prop- 
erties that turbulence must have to transport angular momentum outwards it is instructive 
to take a step back from the MRI and rethink the basic meaning of accretion disk equations. 
Any axisymmetric model that includes turbulent transport is implicitly a mean field theory 
in that locally, turbulence breaks axisymmetry. Therefore equations that evolve quantities 
presumed to be only a function of radius R are only relevant if they can be derived from a 
plausible (even if crude) mean field formalism. 

Here we revisit the hydrodynamic accretion disk equations from a mean field point 
of view and evaluate whether the standard presence of turbulent transport as merely an 
enhanced viscosity emerges from plausible assumptions from a mean field theory. We find 
that it does not. Instead we find transport terms in the surface equation and in the angular 
momentum transport equation that are missing from the standard a viscosity formalism. 
This leads to conditions on the properties of these transport terms that allow outward 
transport of angular momentum. That only turbulence with specific properties can transport 
angular momentum outward is suggested by the simple fact that specific angular momentum 
in a Keplerian disk increases with radius and therefore one might initially expect turbulent 
diffusion to transport angular momentum inwards. In at least some systems this expectation 



is confirmed, as studies of convective turbulence show this effect ( IRyu fc Goodman! Il992 
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Klev et al.lll993f) and ind eed one of the virtues of MRI is that it is expected to drive accretion 



( IBalbus fc Hawleylll998f ). 

In Sec. 2 we derive our mean field formalism and derive the evolution equations for the 
mean surface density and mean orbital angular momentum. We compare the results to the 
conventional versions of these equations in Sec. 3. In Sec. 4 we discuss the implications 
of the new equations for outward angular momentum transport and steady state disks. We 
conclude in Sec. 5. 



2. Mean Field Formalism 

We distinguish three relevant time scales for the accretion disk: (1) a long global time- 
scale t g corresponding to accretion of the entire disk; (2) a dynamical time-scale td < t g 
corresponding to the local accretion time-scale; (3) a short time-scale r <C td corresponding 
to the energy containing eddy turnover time. We assume that that the disk is thin (height 
h <C R) and that the time averaged disk is equivalent to an axially symmetric and vertically 
integrated disk. While full vertic al averaging is appropriate for extracting fluxes, some 



physics is lost (e.g. lBlackmanll200ll ) . Accordingly, we adopt cylindrical coordinates (i?, 0, z). 
We assume that the largest turbulent eddies are small compared to the radial scale at all 
radii, that is, A C i?, where A is the eddy scale. We then decompose all quantities into 
fluctuating and mean values such that q = q + q where q is a flow quantity where q = (q) 
represents the mean and q is the fluctuation about the mean. The fluctuations average to 
zero ove r time-scales l onger than r and the means have temporal variation scales much larger 



than r flMoffattlll978r ). 



In addition, we make the important simplifying anzatz that fluctuating quantities can 
be separated into a factor that varies on spatial scale R times an isotropic homogeneous 
fluctuation. For the surface density and the velocity in particular, we write 

£ = and v = v*u. v (1) 

where v* = (v 2 ) 1 / 2 , u„ is a statistically isotropic and homogeneous vector, = (E 2 ) 1 / 2 
and «2 is a statistically homogeneous scalar. This decomposition is similar to a WKB 
approximation and will prove useful in section 3. 

Using the above decomposition, we now produce a mean-field derivation of the evolution 
of a conserved quantity with areal density q, the relevant examples of which are surface 
density E or orbital angular momentum density L z . For such a conserved quantity, we have 

d -> d _ d _ -> 

— q + V • (qv) = = — q + — q + V • (qv + qw + qv + qv). (2) 
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Taking the time-average (as described above) of Eq. ([2]) we find 

^g + V-(gv) + V-(gv) = 0, (3) 

where transport terms akin to diffusion or viscosity derive from the term V-(gv). Subtracting 
Eq. © from Eq. © we find: 

d - _ 

_ g = -V- (gv + gv + gv- (gv)). (4) 

We approximate the implications of Eq. (jl]) at a point p by estimating the time integrated 
fluctuating flux of q through the surface S composed of those points a distance A = v±t from 
the point p and approximating it as constant for a time r. The distance A is determined at 
the surface S and is chosen because it encloses the volume V that the point p can exchange 
material with over a time r. This allows us to write Eq. (j3J) as: 

q — -— / (gv + gv + gv - (gv)) -dS ~ — — V • A n_1 (gv + g v + q v - (gv)). (5) 
V Js" v ' V 

In the above Eq. (jSJ) F 9 is the fluctuating flux of g as given in Eq. and n is the number of 
spatial dimensions, with n = 2 for our thin disk. The factor A"^ 1 in the right-most term in 
Eq. dSJ) comes from the fact that under our definitions, S oc A n_1 . In producing Eq (jlj), our 
use of a correlation time r has therefore lead us to construct an effective divergence given 
by V-fl U c = A"^ 1 ) V • A"" 1 that is centered on p not in units of distance, but rather in units 
of correlation lengths. Note that this generalized divergence reduces to the basic divergence 
only if A is independent of position. In the case of an accretion disk, A does depend on 
position and varies on the scale R. Upon applying n = 2 to Eq. (jSJ) and using A = v*t we 
then obtain: 

T -» ._ „ 1- _ 

q ~ — - V • A(gv + gv + gv — (gv)) = V • Tv*(q\ + gv + gv — (gv)). (6) 

A ' v* 

We can use Eq. (EJ) to find 

V • (gv) = -V • (u„ ^V^ [tu*(^^+ gv + Wp}) (7) 

inner V- (i) (ii) (in) 

where we have not written the term V • (u„V-(gv)) which vanishes as it is the average of 
fluctuation times a mean. We now discuss how (J7J) simplifies for our quantities of interest: 
surface mass density £ and angular momentum surface density L = L z = Y,R 2 Q. 

For the case q = S, we first note that S„ C S because turbulent density fluctuations 
do not dominate for a non-self gravitating disk with r <C id, and accretion disks have inner 
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boundaries of no return. More specifically, in a near-steady state for which td ^t g , Eq. ([3]) 
implies that 

V- (Ev) ~ V- (Ev). (8) 

Only the radial derivatives contribute so only the radial mean velocity contributes to the 
divergence, and thus Vr/v* ~ £*/£ <C 1. This implies term (i) « term (ii) and also 
that both can be dropped, as they are small with respect to term (Hi) in ([7]) when q — E. 
Then, upon using the chain rule for the inner derivative on term (Hi) the contribution 
that comes from the derivative operating on the fluctuating v vanishes because it involves 
a vector correlation of fluctuating quantities (products of u vectors) which vanish in our 
decomposition. Only the contributions to term (Hi) with the derivative operating on E 
contributes. Finally then, we can write: 

V • (Ev) ~ -V • (Vtv 2 E) ~ -V 2 (z/iS), (9) 

where v\ ~ rv 2 is the diffusion coefficient. The diffusion coefficient v\ is due to the correla- 
tions of fluctuating values. 

We now consider the case q = L = £i? 2 f2 in (J7j). First we note that term (i) can again 
be neglected with respect to term (ii) as this does not depend on the choice of q. To foster 
evaluation of the contributions from (ii) and (Hi) we note that the mean and fluctuating 
values of L are given by 

L = (mR 2 ) = WlR 2 + (£%R) ~ WlR 2 (10) 

and 




respectively. Note that the fourth term in ( flOi) is smaller than the third term for the same 
reason that in (fTTl) term (c) (a) and term (b) (a), namely because v -C Rfl (thin disk) 
and E <C E as noted in the previous paragraph. Then, by analogy to Eq. (P), term (Hi) of 
Eq. <£7§ becomes 

-V 2 (z/il). (12) 

Because the transport coefficient deriving from term (Hi) of Eq. (j7|) is independent of the 
transported mean quantity q, it is plausible to conclude that the transport coefficients of 
Eqs © and (Tl2|) are the same, explaining why we have used v\ in both equations. 

Unlike the case of surface density, term (ii) of Eq. (jTj) can be significant for angular 
momentum. Inserting the fluctuating angular momentum density (Eq. ITTj) into term (ii) of 
Eq. (ED we find 
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The dependence of term (b) of Eq. (fT3]) on v$ provides a preferred direction, so the averaging 
does not eliminate the term in which the inner divergence operates on in (b). In fact this is 
the dominant contribution from (b) since fluctuating quantities vary on scale A <C R. Using 

— * - 

V ~ — , we can write this contribution from term (b) as 

V- (uy^TvlEv^R— >, (14) 

where we have defined /i = v^/v* which parameterizes the dependence on the angular dis- 
tribution of the velocity fluctuations. The significance of the parameter /x can be seen by 
contemplating the difference between radial and azimuthal turbulent forcing: in a Keplerian 
disk radial forcing will move material to positions where it has below Keplerian angular 
momentum while azimuthal forcing, combined with orbital mechanics, will move material to 
positions where it has above Keplerian angular momentum. The averaged quantity in Eq. 
(JHJ), like all mean values, depends only on R so only the radial component of u v survives. 
However, \i depends on v^. Since orbital motion converts into a radial velocity on a time- 
scale Q~ l , only a fraction ~ rf2R of contributes to the correlation. Note that Qt is the 
inverse Rossby number for the disk. Using this in (|14l) . term (b) of Eq. (|13|) can therefore 
be written _ 

r, 1 _™^\ ^ ( "2 L i 



V • UrviZv^R rfiR) = V ■ -=-R , (15) 

V*T ' \ R J 

where u 2 ~ tv\ and includes the system's dependence on the angular distribution of the 
velocity fluctuations (through the dependence of v 2 on //, which is not linear because of the 
averaging process). The quantity z/ 2 can become negative if the velocity fluctuations along 
the orbital motion correlate sufficiently strongly with inwards fluctuating motion. Further, 
while we expect |^i| ~ |i/ 2 |, we do not expect V\ = v 2 . In §3 we will explore the consequences 
of different relative values of u% and v 2 . 

We have yet to evaluate term (a) of fflBl . but we now argue that it is much smaller than 
term (6), the latter contributing Eq. ffl~5l) to our eventual conservation equations. We need 
only compare V • Si? 2 f2v, (a) to V • T.Rv^v, (b). We note from (CQ) that the divergence on 
(a) contributes only a surviving term that pulls out a factor of ~ l/R. Further, we estimate 
that v ~ f ~ v*±. We can then use Eq. (jSJ) to find E ~ Tw/v^ ~ EX/R. We find for (a): 

£R 2 Q ZvRtt SAri?fi _ „ 

~ ~ ~ SAii. (16) 

R U* it 

To estimate (6) we note that the inner divergence pulls out a factor — where the top comes 
from the azimuthal to radial velocity conversion factor and the denominator comes from the 
scale over which an azimuthal turbulent speed varies; note here that this term results from 
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Vcj, which picks out a preferred direction violating isotropy and keeping the spatial scale of 
variation v+t rather than R as in term (a). We thus have for term (b) 

^ Til ~M (17) 

Since X R, (I17p dominates (TIB]) so that term (6) dominates term (a) of Eq. (1131) . 

We can now finally combine Eqs ffl2|) and ffl5|) to find that Eq. ([7]) with q = L gives 

V • (Lv> ~ -V>!l) + V • f^R-J , (18) 
where only radial derivatives are non- vanishing. 



3. Conservation Equations and Comparison to Conventional Form 

We can use Eqs (Q and ffl8|) to write the full, mean-field conservation equations for 
mean surface density £ and orbital angular momentum surface density L = HR 2 Q: 

The mass flow can be characterized as a turbulent flux (right side of Eq. [T9l) and a fall-back 
flux due to orbital angular momentum mismatches (vr). For Q cx R~ 3 ^ 2 and dVt/dt = 0, we 
eliminate vr by combining the above two equations to find the net outwards mass flux 

F M , net = J^E) - v R T, = ±(y&) + ^ - ^^2~L). (21) 

A finite ^ > is important in reducing angular momentum mismatches between turbulently 
transported material and the local material at the new position. 

Eqs. f fT§ll2"Tj) can be contrasted to the analogous equations presented in standard for- 



mulations for Q cx R~ 3/2 flFrank et all 120021 1: 



l E + 14 ( ^ E) = (22 > 



(23) 
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and _ 

- - W (24) 

where i/ is the Shakura-Sunyaev turbulent viscosity. These equations differ from (I19H21I) in 
that there is no turbulent transport term on the right of (Tl9|) and there is only one transport 
coefficient v in the right of (j23|) as compared to the two transport coefficients in ([20]) . The 
transport term in the surface density term (Tl9|) is particularly noteworthy. The presence of 
this term highlights that the mean radial velocity is not the only contributor to mass motion; 
there is also a turbulent diffusion of mass. This is general expected in a turbulent flow and 
so its absence in the standard formalism raises concern. 



What is the origin of the differences between (11911211) and (I22H24P just presented? The 
usual derivation of the conservation equations does not naturally follow from a mean field 
derivation as we have shown by our derivation of (U9H2ip . In particular, the standard ap- 
proach to obtain (!23j) is to simply replace the micro-physical viscosity in the Navier-Stokes 
equation with a turbulent viscosity, assume axisymmetry and vertical integration and then 
derive the angular momentum conservation equation from the momentum density equation. 
Doing so does not produce any transport term in the surface density equation. A key point 
is that although we have written over-lined quantites £ and L in ( 12211241) we do so only 
because it is only the mean quantities that have the axisymmetry and vertical independence 
usually assumed; formally a purely radial dependence requires a mean field formalism We 
have therefore have identified an inconsistency with the standard accretion disk formalism. 

Another expression of the inconsistency is to suppose that a replacement of micro- 
physical viscosity with a turbulent viscosity were a legitimate closure such that no other 
transport terms appeared and the equations remained the same. This would imply that the 
mean field Navier-Stokes equation and turbulent viscosity could be derived via integration 
of suitable Boltzmann equation by analogy to the derivation of the standard Navier-Stokes 
equation and micro-physical viscosity. However, the derivation of the usual Navier-Stokes 
equations comes from the truncation of an expansion in the mean free path divided by the 
macroscopic gradient scale, or a collision time divided by a macroscopic evolution time. 
For accretion disk turbulence, the latter ratio would be replaced by tQ, a quantity not 
guaranteed to be small, and frequently assumed to be of order unity. This highlights why 
a simple replacement of the micro-physical viscosity by a turbulent viscosity is at best, 
incomplete. 
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4. Conditions for Outward Angular Momentum Transport and Steady-State 



Here we investigate some consequences of (U9H2ip . The equation for the net mass flux, 
Eq. (12T1) . can be rewritten as 

F ™»>> = Tnim^"^ - ^wm^). 

This form allows us to explore the consequences of different transport coefficients v\ and v 2 
and different radial density profiles. 

We use the standard formula for the disk density scale height h = c s /Q, where c s is 
the sound sp eed. We scale our transport coefficients to c s h as in a Shakura-Sunyaev a disk 



prescription (IShakura fc Syunyaevl Il973l ). ~ X t v t = aic s h, where the i refers to 1,2 of 
the previous section. In a flared disk whose temperature profile is determined by stellar 
heating (T oc R^ 1 ^ 2 ), we have c s oc i?" 1 / 4 and h oc i? 5//4 . If all «j are constant, then z/j oc R. 
From Eq. 0251) we see that £ oc R~ 3 ^ 2 then results in no net mass transport (although the 
turbulent and fall-back (v R ) terms of e.g. (1211) will in general be separately non-zero), and 
deviations from log^S = —3/2 will result in differently signed mass fluxes (log^ being the 
base R logarithm). For a non-zero steady state mass flux with ui oc u 2 oc c s h, Eq. fl25l) 
would imply log i? (^S) = . In that case, the mass flow (positive being outward) becomes 

M = 2irRF M) net = 2tt - v 2 y)j = 2tt (y - z/ 2 ) S, (26) 

and we need v 2 > V\j1 > to generate an inwards mass-flux. Finally, writing Vi oc fjAt, At = 
v t r and requiring that v t oc c s the condition log^z/E) = becomes log fi (Srfi) = \og R (Q / c 2 s ) . 
For the aforementioned flared disk, log R (Q/c 2 ) = —1. 

We find another implication by noting that v 2 7^ implies the presence of orbital 
oscillations from fl20l) . As long as the amplitudes are small, the oscillations can be treated 
as simple harmonic oscillators and we therefore expect oscillating quantities such as S(i?) 
to be proportional to sin(rfi). For S(i?) to vary as a power law in R, either a constant 
inverse Rossby number tQ (eddy life time scaling with orbit time) or a small rfl such that 
sin(rfi) ~ tQ is then required. A short eddy lifetime is interesting as it lowers v for a given 
v t . (We note that tQ > 1 need not be considered because any eddies with this property 
initially will be sheared such that such that r — > Q^ 1 ; shear makes an otherwise long lived 
eddy shred on a rotation time). 
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5. Conclusion 

Any turbulent accretion disk model in which axisymmetry in surface density is assumed 
is unavoidably a mean field model and therefore the equations governing such a model should 
formally be derived from a mean field theory with a plausible closure. We have presented 
an explicit mean field approach for deriving the mean surface density and orbital angular 
momentum transport in accretion disks and found that two transport coefficients rather than 
just one emerge most naturally. Notably we find that the conventional a formalism misses 
a turbulent diffusion term in the surface density equation. 

In constraining the new transport coefficients for specific cases of inward, outward, and 
zero net mass transport, we find that: (1) log i? (z/S) = and v 2 ^ ui/2 is the condition for a 
steady-state disk with non-zero net mass transport with v 2 > vi/2 the condition for inwards 
net mass transport and (2) \og R {yll) = —1/2 or v 2 = V\j2 are the conditions for a steady- 
state disk with zero net mass-transport; (3) For a finite is 2 , the inverse Rossby number tQ 
must be either constant or small for disks to have monotonic power-law radial dependences of 
E. This emerges because a finite v 2 implies oscillations in the angular momentum equation. 
Such a condition on rQ is not unexpected if one considers that these oscillations depend on 
sin(rf2) and since rfl ~ 1 blurs the distinction between fluctuating and orbital time scales, 
significant radial variations in tQ would have implications for mean surface densities and 
fluxes. 
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